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Abstract 

We introduce a model that simulates a kinetic roughening process with two 
kinds of particles: one follows the ballistic deposition (BD) kinetic and, the 
other, the restricted solid-on-solid (KK) kinetic. Both of these kinetics are 
in the universality class of the nonlinear KPZ equation, but the BD kinetic 
has a positive nonlinear constant while the KK kinetic has a negative one. In 
our model, called BD-KK model, we assign the probabilities p and (1 — p) to 
the KK and BD kinetics, respectively. For a specific value of p, the system 
behaves as a quasi linear model and the up-down symmetry is recuperated. 
We show that nonlinearities of odd-order are relevant in these low nonlinear 
limit. 
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The growth of interfaces by nonequihbrium kinetic roughening models is a very inter- 
esting topic of the far from equihbrium statistical mechanics [1-3]. In the two past decades, 
several models have been proposed, as ballistic deposition |Q, Eden model 0, solid-on-solid 
(SOS) model with surface relaxation and SOS model with refuse SOS model with 
diffusion 0. In computer simulations, interfaces are described by a discrete set {hi(t)} 
which represents the height of a site i at the time t. Such interface has L'^ sites, where L is 
the linear size and d is the dimensionality of the substrate. The roughness of the interface 
is defined by 

u;\L,t) = /^-^^^ih,-hr'^ , (1) 

where h is the mean height at the time t and the symbols < .. > means average over inde- 
pendent computational samples. In most of the kinetic roughening models, the roughness 
obeys the Family- Vicsek dynamic scaling 



co{L,t)^L-f(^-^^ , (2) 

where the function /(x) must be L-independent. The roughness behaves as ~ t^, for short 
times (1 ^ t ^ L^) and behaves as uJoo{L) ~ in the steady state. The f3 and a are the 
growth and roughness exponents, respectively, and are related with the dynamical exponent 
z through the relation z = a/ (3. For some systems, a = /5 = and z 7^ which means that 
the roughness does not obey Eq.2 and has a logarithmic behavior in space and time. 

The kinetic growth models are also described by the continuum Langevin-like equations in 
the coarse-grained limit. These equations have terms which represents the main interactions 
among the incoming particles. Example of a linear equation is the EW equation 0, 

(9/i(x,t) 



dt 



Vo + v{^,t) + iyV^h{x,t) , (3) 



that describes the fluctuations of the SOS model with surface relaxation |^. In Eq.3, the 
first two terms of the right side are related to the deposition of particles. This deposition 
has a rate vo and a noise with zero mean and variance given by 



ri{x,t)ri{x,t)) = D5'^{x-^)5{t-t) . (4) 

The third term represents the surface relaxation process. 

The exponents of Eq.3, obtained by Fourier analysis |§|,|TT|, are P{d) = (2 — (i)/4, a{d) = 
(2 — d)/2 and z{d) = 2. For d = 1, these expressions gives (3 = 1/4 and a = 1/2. For d = 2, 
the scaling exponents are P = a = 0. 

There are also some nonlinear equations that describes nonlinear kinetic roughening 
models. The more known nonlinear equation is the KPZ equation [|l^] 

^^^ = z;o + r;(x,t) + z/V^Mx,t) + ^(VMx,t)f . (5) 

This equation is more complete than Eq.4 because the nonlinear term may represent the 
lateral growth, or the appearance of a driven force. In ci = 1, the exponents of this equation 
|l^ are f3 = 1/3, a = 1/2 and z = 3/2. In = 2, the analytical solution is not known. 



Examples of models in the universality class of KPZ equation are the ballistic deposition 
(BD) and the SOS model with refuse (KK) |^. In d = 1, numerical simulations |T^ 



indicate p ^ 0.30, a ^ 0.47 for the BD model and p ^ 0.332, a ^ 0.489 for the KK model. 
In ci = 2, these exponents are P ^ 0.24, a ^ 0.40 (BD) and p ^ 0.25, a ^ 0.40 (KK). 

In this letter we report on results of computer simulation of a growth model with two 
kinds of particles. Both of them obeys the KPZ kinetic, bu they have opposite signs of the 
nonlinear A constant. The effort to understand the nonlinearity in stochastic systems out 
of equilibrium is due to the great influence of the nonlinearities in the scaling analisys, as 



pointed by Binder et al [14|. The anisotropic KPZ equation, studied analytically in 1991 by 



Wolf [T^, also contains two nonlinear terms with opposite signs that describes the distin- 



guishable directions of the substrate in a vicinal surface. He has found, by renormalization 
group, a logarithmic behavior of the roughness {a = P = Q) for opposite signs of nonlin- 
earity and power-law behavior when the nonlinearities have same signs. This results were 
confirmed by Kim et al. |T6| by simulations of a discrete model with two different kinetics 
applied in each direction of the substrate. 



3 



The motivation of our work is the paper pubhshed by Bernardes et al [0. They have 



studied the deposition of particles with different radii on a cold substrate by Monte Carlo 
simulation. The authors have found the growth exponent (3 ~ 0.26, for d = 1. So, they have 
concluded that the universality class of this model is next to the EW class 0. However, 
there are two nonlinear characteristics in the morphology of the model of Bernardes et al: 
(i) existence of porosity in the bulk; (ii) the growth velocity, that is, v = (dh/dt^, is greater 
than the deposition rate. These characteristics can indicate that the up-down symmetry 
{h — s> —h) is broken The break of this symmetry leads to appearance of the nonlinear 
term in the KPZ equation. Bernardes et al. also have checked the need of logarithmic 



correction in the behavior of roughness ||l^ , p!8| that indicates the presence of odd nonlinear 
terms in the growth equation. 

The aim to create a model with opposite KPZ nonlinearities is to verify the possibility 
of generating a linear process (when the effective nonlinearity vanishes) with a morphology 
of a nonlinear growth process. Our model is a probabilistic combination among the ballistic 
deposition model (BD) and the SOS model with refuse (KK). The KK model occurs with 
p-probability, and the BD model occurs with (1 — p) -probability. 

Ballistic deposition is a process where particles are dropped vertically onto a smooth 
substrate. The incoming particles are automatically joined to the growing cluster when its 
first contact with the growing interface occurs. In the in-lattice version, we select at random 



a site i of the lattice and its new height is evaluated by the algorithm [19 



h[ = max{hi + 1; , (6) 

where {j} are the first neighbors of the site i. The BD kinetic does not generate a solid- 
on-solid deposition because it generates a structure with porosity. Therefore, we define the 
growing profile as the major height of an occupied site of each column. 

The KK model is a SOS random deposition with the difference of height constraint 
hi — h^jj < m, where m is the parameter that controls the roughness. If the height of the 
particle deposited on the site i did not satisfy the height constraint, this particle is not 



incorporated to the interface. 

Both of them are in the same universahty class of KPZ equation. We have chosen these 
two models because: (i) the BD model generates a bulk with porosity and has nonlinear 
parameter Abd > because the growth velocity is bigger than the deposition rate vq] (ii) 
the KK model has Xkk < because the kinetic of refuse makes the growth velocity smaller 
than the rate of deposition. 

In our simulations, a unit of time means that we have done L attempts of deposition. 
Moreover, all simulations were done with a one- dimensional substrate (d=l) and, in the KK 
model, the difference of height constraint m = 1. 

Figure 1 shows the plot of the effective growth exponent /3e// vs. the parameter p, for 
L = 50, 000. The long dashed line is the exact value of (3 obtained from KPZ equation by 
renormalization group. We have done 100 independent runs for each probability and we 
applied consecutive slopes method |0] in the log-log plots of u vs. t in the time interval 
20 < t < 10000, giving an ensemble of jSeff exponents for each value of p. So, we have 
estimated the error bars around each value of Peff- For p* = 0.83, the model is near to the 
EW class, because the growth exponent is /?e// ~ 0.27. However, the error bars in this region 
has increased, indicating the need of more careful analysis of the scaling. These results can 
point that the effective KPZ term was removed in the probability p*. 

In order to make better characterization of the KPZ nonlinearity in our model, we do 
finite size analysis in the growth velocity. In 1990, Krug and Meakin pO[ showed that the 



finite size correction, Av{L,t) = v{L,t) — v^o, for a model in KPZ class behaves as 

Av{L) ~ -AL-"ii , for t > , (7) 

where the a\\ exponent depends on the roughness exponent. The Av correction goes to zero 
when the KPZ term vanishes. So, with Eq.7, we can obtain the sign of the KPZ nonlinearity 
and determine when the nonlinearity goes to zero in function of the tunning parameter p. 
Figure 2 shows the plot of Av = v{L = 10) — v{L = 1280) vs. p for the BD-KK model. The 
finite size correction vanishes for p ^ 0.81 (see inset), very close to the value obtained in the 



minimum of the (5 vs. p plot (see Figure 1). 

We also check the need of multiplicative logarithmic corrections in the scaling for p = p*. 
When the growth equation for a process has a sequence of odd nonlinear terms as 



for n = 1, 2, the roughness behaves as jl^JTsU 



co{L,t)r^t'/\\ogtY/' , (9) 

for t <^ U . So, if logarithmic corrections are accepted for p = p*, it means that the system 
is marginally in the EW class. 

In order to show that Eq.9 is really the better equation that describes the system near 
to p*, we rewrite this equation as 

u;{L,t)r^t\\ogty , (10) 

and we do small variations around the exact values 5 = 1/4 and 7 = 1/8. We analyze 
the validation of the Eq.lO by the evaluation of the deviations from the horizontal curve 
Y{6,-f,t) = uj{L,t)/t\logty' vs. t, because this emphasizes better the deviations of the 
behavior from this equation. So, we measure the relative error AY/ {¥) of each curve 
Y{6,'j,t) in function of the variations in 6 and 7. Figure 3a shows the plots of Y{t,'~f) vs. t 
with 6 = 1/4 and Figure 3b shows the plots of Y{t, 6) vs. t with 7 = 1/8. The insets shows 
the relative error AY/ (Y) vs. the scaling exponent related to each case, 7 or S, respectively. 
The relative error, for the two cases, reaches a minimum when the exponents 7 = 1/8 and 
6 = 1/4, indicating that Eq.9 is a good description for the temporal behavior of roughness. 
This is pointing out that odd- nonlinear terms are relevant. 

The nonlinear Eq.8 preserves the up-down symmetry, but this is not obvious in BD-KK 
model at p = p*, and also in the growth model of Bernardes et ai. Eq.l can be generalized 
for any moment q of the height distribution, in d = 1, as 

^''(^'^) = (^E(^^-^)') ' (11) 



and we concentrate on the behavior of odd-moments, especially the third moment {q = 3) 
which is related to up-down symmetry. The skewness, defined hy S = cu^ / (co'^y^'^ , can 
show us if the system has this symmetry. If the skewness is null, the system has up-down 
symmetry because all odd moments of the height distribution vanish. On the other hand, for 
systems without the up-down symmetry and in the KPZ class, some authors [PT],p2| believe 
that the skewness has a universal value |iS| ~ 0.28. Figure 4 shows the skewness iS as a 
function of the time t for the model with positive nonlinearity {p = 0.0, triangle up), with 
negative nonlinearity {p = 1.0, triangle down) and at the low nonlinear point p = p* (filled 
circles). As illustration, we also show the curve for the model with surface relaxation 
which is a linear model and has the up-down symmetry (plus). We note that the skewness 
for p = tends to iS = 0.28 and, for p = 1.0, to S = —0.28. For p = p*, in the asymptotic 
limit, the skewness goes to zero, suggesting the presence of the up-down symmetry. 

In conclusion, we have studied the scaling properties of a model with opposite signs of 
the KPZ nonlinearity through numerical simulations of a model with two kinds of particles. 
The Kim-Kosterlitz kinetic occurs with a probability p and the ballistic deposition model 
with a probability (1 — p). For a specific value of the tunning parameter p, we show that 
the KPZ nonlinearity goes to zero and the up-down symmetry is recuperated. We also show 
that odd-nonlinearities are relevant in this model. 
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Americo T. Bernardes for fruitful discussions. The numerical simulations were made in part 
in Alphas 500 Au of the Departamento de Ffsica -UFMG and in part in a Sun HPC 10000 
of the Cenapad MG-CO. This work was supported by CNPq, Fapemig and Finep/Pronex, 
Brazilian agencies. 
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Figure captions 
Figure 1 

The growth exponent /3 vs. the parameter p plot of the BD-KK model for L — 50, 000. 
The long-dashed line is the exact value of the /3 exponent for the KPZ equation. 
Figure 2 

Plot of the difference between the steady state growth velocities for L = 10 and L = 1280 
vs. the parameter p which gives the amount of KPZ nonlinearity in the system. The inset 
is showing the behavior next to the crossover. 

Figure 3 

Plots of the function Y{t,j,d) vs. the time t: (a) 6 = 1/4 and small variations in 7 are 
performed, (b) 7 = 1/8 and small variations in S are performed. The bold curves are showing 
the behavior of the function Y(t, 1/8, 1/4). The insets shows the relative error AY/ (Y) vs. 
the scaling exponent related to each case. 

Figure 4 

Plots of the skewness 5 as a function of the time t for p — 0.0 (triangle up), p — 
1.0 (triangle down) and p = p* (filled circles). The skewness for the EW linear model is 
represented by the symbol plus. The long dashed straight hues indicate the ±0.28 estimated 
values. All simulations were done with L — 50, 000. 
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